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Abstract 

Comment on "Classification Scheme for Plienomenological Uni- 
versalities in Growth Problems in Physics and Other Sciences" by P. 
Castorina, P. P. Delsanto and C. Guiot, Phys. Rev. Lett. 96, 188701 
(2006) is presented. It has been proved that the West-like function 
of growth derived by the authors is incorrect and the approach does 
not take into account the growth of the biological systems undergoing 
atrophy or demographic and economic systems undergoing involution 
or regression. A simple extension of the model, which permits deriva- 
tion of the so far unknown involuted Gompertz function of growth is 
proposed. 

The concept of phenomenological universalities introduced recently by 
Castorina, Delsanto and Guiot (CDG) [1] is a useful tool for investigation of 
the nonlinear processes in complex systems whose dynamics is governed by 
the system of equations 

f?^ = .M,M ^ = -tW. (1) 
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Here, r = x{0)t denotes dimensionless temporal variable, whereas is 
a generating function, which extended into a series of x-variable generates 
different functions of growth for a variety of patterns emerging in complex 
systems in physics, biology and beyond. For example for $ = a; one gets 
the Gompertz function, whereas for $ = a; + bx'^ the allometric West-like 
function is derived [T] 



yiT)G = exp [1 - exp(-r)] , y{T)w = [l + b- bexpi-r)]'^^'-'^ . (2) 



Unfortunately, the West-hke function obtained by CDG [T] is incorrect as it 
is not a solution of the differential equation (9) in [1]. Additionaly in the limit 
limt^o 1/(t)vk = 1 it does not produce the Gompertz function as indicated 
by the authors [JJ. To explain those inconsistencies Eqs. ([T]) were solved for 
$ = X + bx"^ employing Maple vs. 7.0 processor for symbolic calculations. 
The calculations provided the correct solution y{T)w = [1 + b — 6exp(— r)]^^'' 
with power 1/b and not 1/(1—6) as derived by CDG [T]. Employing the Maple 
one may also prove that 



as it should be. 

Although CDG claimed [1] that they have developed a simple scheme 

that allows the classification of all the growth problems described by Eqs. 

this approach does not take into account the growth of the biological 
systems undergoing atrophy or demographic and economic systems undergo- 
ing involution or regression. In biological systems such a situation appears 
in avian primary lymphoid organs: thymus and bursa of Fabricius as well 
as in thymus of mammalians. To extend the CDG approach and derive an 
involuted function of growth, we employ a linear expansion of the generat- 
ing function = CiX + Cq, which includes a constant term cq omitted 
in the CDG scheme and Ci coefficient, which in the previous approach was 
constrained ci = 1 . Additionally, we assume that x{t = 0) = (1 — co)/ci, 
which for cq = and Ci = 1 gives the CDG condition x(0) = 1. Employing 
the relationships ([1]) one gets the involuted Gompertz function (cq, ci > 0) 



which can be specified in the form applicable to direct fitting of the experi- 
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mental data 



y{t) = yoexp |- [1 - exp(-at)] | exp {-bet) . (5) 

Here, we use the correspondences cf = a/b, Cit = at and c = Cq. To prove 
that function (jS]) correctly describes the evolution and involution of organs 
undergoing atrophy, we employed it to fit 30 mean absolute weights of the 
thymus of male Wistar rats evaluated in the period of 1-780 days. In the 
calculations we used thymuses sampled from strain rats of the mean age of 
24 months, originating from randomly mated culture (Department of Toxi- 
cology, University of Medical Science in Poznah, Poland). The mean values 
of the absolute thymus weight were fitted to the 4-parameteric function ([5]) 
using the weighted nonlinear least-square routine with statistical weights 
taken as inverse squares of uncertainties u{i) of the thymus weights. As the 
criteria of goodness of the fit, the standard deviation a and normalized stan- 
dard deviation a were employed. The calculations provided the values of 
the parameters yo = 0.0105(16) [g], a = 0.0472(44) [1/day], b = 0.1914(194) 
[1/day] and c = 0.0140(38) [1/day] reproducing the experimental data with 
a = 0.011 [g] and a = 1.64. All parameters fitted are statistically well evalu- 
ated and the correlation coefficients between them cover a satisfactory range 
(0.6455 -0.8971). 

The results obtained indicate that the Ui class introduced by CDG [T] 
should be divided into two subclasses Un - representing the sigmoidal Gom- 
pertzian growth, and f/12 - including the involuted Gompertzian growth. 
The f/12 solution ([5]) is a very important contribution to the classes obtained 
previously as it is usefull to fit the data for the systems undergoing atrophy. 
The sigmoidal (S-shaped) functions (Gompertz, West, logistic, Richards etc.) 
cannot be apply to this aim as they correctly describe the ideal situation in 
which exponential growth is exponentially retarded and saturated as time 
continues. 
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